NeXSPheRIO results on azimuthal anisotropy in Au-Au collisions at 200A GeV 



Wei-Liang Qian,^ Rone Andrade/ Arthur dos Reis/ Yogiro Hama,^ Frederique Grassi,^ and Takeshi Kodama^ 

^Instituto de Fisica, USP, C. P. 66318, 05314-970 Sao Paulo, SP, Brazil 
^Instituto de Fisica, UFRJ, C. P. 68528, 21945-970 Rio de Janeiro-RJ , Brazil 

In this work, we present the results obtained by the hydrodynamic code NeXSPheRIO on 
anisotropic flows. In our calculation, we made use of event-by-event fluctuating initial conditions, 
and chemical freeze-out was explicitly implemented. We studied directed flow, elliptic flow and forth 
harmonic coefflcient for various hadrons at different centrality windows for Au-|-Au collisions at 200 
AGeV. The results are discussed and compared with experimental data from RHIC. 
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PACS numbers: 25.75.-q, 24.10.Nz, 25.75.Ld 



I. INTRODUCTION 

Anisotropic flows are defined by the Fourier decom- 
position of the azimuthal distribution of particles pro- 
duced in heavy ion collisions with respect to the event 
plane. Directed flow and elliptic flow are quantified by 
the first harmonic (ui) and the second harmonic (^2), 
respectively. The study of anisotropic flows in high en- 
ergy nuclear collisions provides us valuable information 
on equation of state (EOS), degree of thermalization and 
the early stage of the created hot, dense matter. It also 
plays an important role in the extensively discussed ques- 
tion whether the quark-gluon plasma (QGP) created in 
Au-|-Au collisions at the Relativistic Heavy Ion Collider 
(RHIC) represents a "perfect" liquid. 

Hydrodynamics offers a reasonable description of the 
experimental data for bulk properties of low pr particles 
for heavy- ion collisions at RHICP, 0, H, 0,Ja,lll- Be- 
sides, the discovery of the "perfect" liquid^, Q at RHIC 
in Brookhaven National Laboratory (BNL) was based on 
the successful modelling of the anisotropic expansion of 
the matter in the early stage of the reaction by means of 
ideal hydrodynamics. Specifically, at the RHIC energy 
of 200 A GeV, the observed V2 values agree well with 
predictions from ideal fluid dynamics [9] in the midrapid- 
ity region {\rj\ < 1) with small or intermediate impact 
parameters and transverse momenta. The agreement be- 
tween model predictions and the experimental data sug- 
gests that hydrodynamic model is well justified to de- 
scribe the intermediate stage of the reaction. It is the 
hydrodynamics that plays an important role to fill the 
gap between the static aspects of QGP properties and 
dynamical aspects of heavy ion collisions measured ex- 
perimentally. 

In this work, we present recent results on azimuthal 
anisotropy obtained by a hydrodynamic code developed 
by our group, called NeXSPheRIO. The code is based on 
full three dimentional ideal hydrodynamics. We calcu- 
lated directed fiow coefficient vi, elliptic flow coefficent 
V2 and forth harmonic coefficient va of Au+Au collisions 
at 200 A GeV at RHIC, for various identified particles 
at different centrality windows. Two important aspects 
are employed in our present approach, namely, fluctuat- 
ing initial conditions (IC) and chemical freeze-out. Their 



effects on observables are addressed. 

In what follows, we will first give a brief description 
of the NeXSPheRIO code and how this code is used to 
compute the observables of our interest. The results are 
shown in the third section, along with discussions and 
perspectives. 



II. THE NEXSPHERIO CODE 

The numerical tool for our present study is called NeX- 
SPheRIO. It is a conjunction of two codes, NeXus and 
SPheRIO. NeXus^' is as an event generator. It is a 
microscopic model based on Gribov-Regge theory [ll|. 
which provides the initial conditions (IC) for a macro- 
scopic treatment of the hydrodynamic stage of the colli- 
sion. The IC state is given as spacial distributions of en- 
ergy momentum tensor and other conserved charge cur- 
rents. We assume, as in most hydrodynamic models on 
the market, that the local equilibrium is reached within 
a short time scale r ~ 1 fm/c. Owing to the partonic 
structure of high momentum incident nuclei, the IC ob- 
tained at a very early instant can not be smooth, but 
highly inhomogeneous in the space. In addition, since 
the size of the system we are dealing with is very small, 
big event-by-event fluctuations are expected in real colli- 
sions. This feature is very different from most hydrody- 
namic treatements where averaged smooth IC are usually 
adopted. As already observed and studied previously in 
[S [H, Q [S m, [13, such fluctuations in IC lead to 
significant effects on experimental observables. In this 
context, one of the advantages of NeXus is that it can pro- 
duce realistic inhomogeneous IC in event-by-event basis, 
therefore provides us with a good opportunity to study 
the effect of fluctuations in IC. We show in Fig.l an ex- 
ample of one such fluctuating event, generated by NeXus, 
for a random central Au-f Au collision at 200AGeV. As 
one can see, the energy-density distribution presents sev- 
eral blobs where energy is highly concentrated. 

To solve the hydrodynamic evolution for systems like 
the one shown in Fig.l is not a trivial task. The section of 
our code which undertakes this job is called SPhcRIO 6j. 
The code is based on the Smoothed Particle Hydrody- 
namic(SPH) algorithm, which was first introduced for 
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astrophysical applications and later adopted in rela- 
tivistic heavy- ion collisions [l9j. The model parameterizes 
the matter flow in terms of discrete Lagrangian coordi- 
nates, called SPH particles. In terms of SPH degrees of 
freedom, the equations of motion can be derived by using 
the variational principle [20|. The main advantage of the 
method is that it is rather robust to deal with any kind 
of geometrical structure and violent dynamics. At the 
meantime, it is reliable and efficient, provided the size of 
SPH particles is appropriately chosen. 

We use in the present study an equation of state (EOS) 
which takes into consideration the strangeness conserva- 
tion. We adopt hadronic resonance model with finite vol- 
ume corrections to describe the matter in the hadronic 
side, where the main part of observed resonances in Par- 
ticle Data Tables [llj has been included. For quark gluon 
plasma phase, the MIT bag model is made use of. By 
solving the Gibbs equation, one obtains a first order 
phase transition in the mixed phase. As a good approxi- 
mation, local strangeness neutrality is assumed. We shall 
neglect any dissipative effects in this work. 

As it was shown [2^. hydrodynamic model without 
chemical freeze-out gives good description of transverse 
momentum spectra of light hadrons. While it reproduces 
the shape of spectra for hyperons and anti-protons, there 
are visible discrepancies in the multiplicities. To com- 
pensate this, we implemented early chemical freeze-out 
for strange hadrons such as A, S, fl and cj) since their 
cross sections are small. At chemical freeze-out tempera- 
ture Tc/i, these particles cease to have inelastic collisions 
and therefore their abundances are fixed. We assume 
complete chemical equilibrium on the surface of chemi- 
cal freeze-out, thus no strangeness saturation factor 75 is 
introduced in our treatment. In our calculation, Tch and 
Tth serve as adjustable parameters as function of cen- 
trality according to experimental data. The introduction 
of chemical freeze-out for hyperons into hydrodynamics 
model helps to correctly reproduce the hyperon yields 
of transverse momentum spectra. In the present work, 
above mentioned chemical freeze-out is employed. 



III. RESULTS AND DISCUSSIONS 

Here, we first tune the free parameters of the model 
to match the experimental data for the multiplicity and 
radial flow. We begin by fixing the IC so as to re- 
produce properly the pseudo-rapidity distributions of 
charged particles in each centrality window. This is done 
by applying an 77-dependent factor ~ 1 16] to the ini- 
tial energy-density distributions for each centrality win- 
dow. Next, to correctly reproduce the transverse mo- 
mentum spectra of charged particles, we appropriately 
parametrize the centrality dependent freeze-out temper- 
ature Tth- As for chemical freeze-out temperatures, they 
are adjusted to reproduce the yield of As. We calculate 



the anisotropic flows of identified particles with respect 
to the event planes of all charged particles. As in our 
previous studies [Tsl [l6l \it\ . the event planes are defined 
using PHOBO subevent method ^23^. 

Our main results are shown in Fig. 2 - Fig. 6. We plot 
in Fig. 2 - Fig. 4 the elliptic flow coefficient V2 of pions, 
protons, anti-protons and As for diferent centrality win- 
dows and minimal bias. In Fig. 5 and Fig. 6, directed flow 
coefficient vi and the fourth harmonic coefficient W4 for 
all charged particles are depicted. The numerical results 
are compared with data from STAR Collaboration 24 1 . 
One finds that V2 and V4 increase with increasing trans- 
verse momentum as well as with more peripheral colli- 
sions, and V2 decreases when the mass of identified par- 
ticle increases. It is observed that the results obtained 
by using sudden Cooper-Frye freeze-out undershoot the 
V2 data of pions and slightly overestimate those of pro- 
tons, anti-protons and As in small transverse momentum 
region with pt < 2, but overall they are in qualitative 
agreement with the experimental data. The result of di- 
rected flow coefficient vi shows reasonable agreement in 
small 7y region. Since STAR data of vi do not extend 
to the regions in which |?7| > 4, we are not able to com- 
pre the data there where the calculated curve turns over. 
Our result of V4 stays little below the experimental data, 
but it demonstrates the correct trend. 

A meaningful solution to improve the elliptic flow 
result is to adopt the continuous emission (CE) 
prescription [26l. \2^. In the CE picture, the emission of 
hadrons occurs not only from the sharply deflned freeze- 
out hypersurface, but continuously from the whole ex- 
panding volume of the system at different temperatures 
and different times. As a result, the large-pT particles 
are mainly emitted at early instants when the fluid is 
hot and mostly from its surface. One might imagine that 
isotropic emission of large pt components at the surface 
will naturally decrease V2 at large transverse momentum 
region, and the px spectra will become more concave. 
These effects are very similar to those discussed in other 
recent studies where one incorporates viscosity explic- 
itly into the equation of motion[28., 29,, 30., 311, O'' con- 
nects the hydrodynamic model to a hadronic transport 
modelfSsj. On the other hand, small-pT components are 
mostly emitted at later times when the fluid is cooler 
and from larger spatial domain, which might lead to an 
enhancement of V2 in small px region. For particles like 
(p, it is usually understood that they have smaller cross 
sections [325, so it is reasonable to treat their freeze-out 
process differently. Consequently, the value of V2 on light 
hadrons such as pions might increase thus the result will 
be improved, and meanwhile those of heavy hyperons will 
remain unchanged. Further work to improve and extend 
the present study is under progress. 
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FIG. 1: Energy density distribution of one random event 
at mid-rapidity plane for central 200 A GeV Au-fAu colli- 
sions generated by NeXus. The diagram is plotted in units of 
GeV/fm^ 




FIG. 2: i;2 of pions as a function of pr at different centrality 
windows; the experimental data are from STAR Collaboration 
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FIG. 3: V2 of protons and anti-protons as a function of pr at 
different centrality windows; the experimental data are from 
STAR Collaboration ■^]. 




FIG. 4: V2 of As as a function of pT at the centrality window 
as indicated; the experimental data are from STAR Collabo- 
ration (24] . 




FIG. 5: of all charged particles as a function of pr; the 
experimental data are from STAR Collaboration [24j . 




FIG. 6: W4 of all charged particles as a function of pr; the 
experimental data are from STAR Collaboration [2^ . 



